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We present recent lattice results on the baryon spectrum, nucleon electromagnetic and axial form 
factors, nucleon to A transition form factors as well as the A electromagnetic form factors. The 
masses of the low lying baryons and the nucleon form factors are calculated using two degenerate 
flavors of twisted mass fermions down to pion mass of about 270 MeV. We compare to the results 
of other collaborations. The nucleon to A transition and A form factors are calculated in a hybrid 
scheme, which uses staggered sea quarks and domain wall valence quarks. The dominant dipole 
nucleon to A transition form factor is also evaluated using dynamical domain wall fermions. The 
transverse density distributions of the A in the infinite momentum frame are extracted using the 
form factors determined from lattice QCD. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 
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I. INTRODUCTION 

During the last five years we have seen tremendous 
progress in dynamical lattice simulations using a num- 
ber of different fermion discretization schemes with quark 
masses reaching closer to the physical pion mass. Many 
collaborations are contributing to this progress. The Eu- 
ropean Twisted Mass Collaboration (ETMC) is using 
twisted mass fermions (TMF), which provide an attrac- 
tive formulation of lattice QCD that allows for automatic 
0{a) improvement, infrared regularization of small eigen- 
values and fast dynamical simulations pj, ^ . Automatic 
0{a) improvement is obtained by tuning only one param- 
eter requiring no further improvements on the operator 
level. A drawback of twisted mass fermions is the 0{a?) 
breaking of isospin symmetry, which is only restored in 
the continuum limit. In the baryon sector it has been 
shown that this isospin breaking is consistent with zero 
within our statistical accuracy by evaluating the mass 
difference between A++(A-) and A+(A°) i3|,i4]. This is 
in agreement with a theoretical analysis [a, Q that shows 
potentially large 0{a?) flavor breaking effects to appear 
in the 7r°-mass but to be suppressed in other quantities. 
A number of collaborations, as for ex amp le QCDSF 0, 
PACS-CS [|, BMW [g^l and CERN fTo| are using im- 
proved Clover fermions for their simulations. It is worth 
mentioning that PACS-CS has simulations very close to 
the physical pion mass albeit in a small volume, whereas 
the Wuppertal group recently calculated meson masses 
and the decay constants using Np = 2-1-1 configurations 
simulated at the physical pion mass [ll|. A number of 
groups adopted a hybrid approach to compute hadronic 
matrix elements taking advantage of the efficient simula- 
tion and availability of staggered sea fermions produced 
by the MILC collaboration 12] and the chiral symmetry 
of domain wall fermions. The Lattice Hadron Physics 
Collaboration (LHPC) has been particularly active in 
producing results on a number of key observables p^lT3|. 
some of which will be discussed in Sections IV and 



V. A very promising recent development is the simula- 
tion of dynamical chiral fermions using large volumes 
and at small enough pions masses. The RBC-UKQCD 
collaboration is generating gauge configurations using 
Np = 2-1-1 domain wall fermions (DWF) [l^, whereas 
the JLQCD Collaboration is producing dynamical config- 
urations with two flavors of overlap fermions [l6l|. Most 
of the current simulations are done using volumes of spa- 
tial length L such that m-^L > 3.5 to keep finite vol- 
ume effects small. The fact that simulations in the chiral 
regime are possible is to a large extend due to algorith- 
mic improvements that yield better scaling behavior as 
the physical pion mass is approached. For a discussion on 
the scaling and a comparison among the different fermion 
discretization schemes see Ref . 1171 . 



II. HADRON SPECTRUM 

The masses of the lowest lying hadrons of a given set 
of quantum numbers are readily calculated by comput- 
ing the two-point function at zero momentum: Ch{t) — 
^,j.(0| J/i(x, t) J^(0)|0). Choosing good interpolating 
fields and applying smearing techniques ensure ground 
state dominance at short time separation t so that gauge 
noise is kept small [1^]. In Figs. [T] and [2] we compare re- 
cent results on the low lying baryon spectrum using dy- 
namical twisted mass 0, Il9j | and clover fermions [J and 
within the hybrid approach (staggered sea and do- 
main wall valence quarks). The level of agreement of lat- 
tice QCD results using a variety of fermion discretization 
schemes seen in Figs. [1] and [2] before taking the contin- 
uum limit or other lattice artifacts into account is quite 
impressive. Small discrepancies seen mainly in the decu- 
plet masses can be attributed to lattice artifacts and a 
systematic analysis of these effects is performed by each 
collaboration before extracting the final continuum val- 
ues. In particular results using staggered fermions may 
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FIG. 1: Comparison of masses for the low lying octet baryons. 
Results using Nf = 2 TMF are shown by the filled (black) 
triangles for L = 2.1 fm and (blue) squares for L = 2.7 fm with 
a — 0.089 fm and with the open (red) triangles for L = 2.1 fm 
and a = 0.070 fm. Results with the hybrid action are shown 
with the (green) asterisks for a = 0.124 fm and results using 
Nf = 2+1 Clover fermions with the open (orange) circles and 
a = 0.0907 fm. For the nucleon we also show results using 
Nf =2+1 staggered fermions (filled (light blue) circles). The 
physical masses are shown by the (purple) star. 



suffer the most from cut-offs effects since the lattice used 
is rather coarse as compared to tliose using twisted mass 
and Clover fermions which have lattice spacings smaller 
than 0.1 fm. 



Having reliable methods to extract the masses of the 
low lying hadrons one can investigate techniques for the 
extraction of the masses of excited states. A number of 
approaches exist. A commonly used method is based on 
the variational approach [2l|: For a given N x N corre- 
lator matrix Cfc„(i) = (0| Jfe(t) 4 (0)|0) one defines the N 
principal correlators with Xk{t,to) as the eigenvalues of 



FIG. 2: Comparison of masses for the low lying decuplet 
baryons. The notation is the same as that of Fig. [T] 

C{to)-^^^C{t)C{to)-'^/'^, where to is small. Since 

lim Afe(<, to) = e-^*-*")^" (1 + e-"^^") , k = l,...,N 

(1) 

the TV principal effective masses tend (plateau) to the N 
lowest-lying stationary-state energies of the hadrons with 
the same quantum numbers. It is crucial to use very good 
operators so noise does not swamp signal and construct 
spatially extended operators using smearing of the quark 
fields as well as applying link variable smearing. The use 
of a large set of appropriately constructed operators is 
also very important. Despite recent calculations using 
this method [2^] the issue of the ordering of the Roper 
resonance as compared to the negative parity partner 
of the nucleon still remains unresolved. Maximum en- 
tropy methods have also be developed for the analysis 
of hadron two-point correlators and recent results can be 
found in Ref . [23] . A new method that relies solely on x^- 
minimization with an unbiased evaluation of errors can 
be applied to extract the masses of the states on which 
the two-point correlator is sensitive on [2^ . This method 
was applied to extract the excited states of the nucleon 
using local correlators that are easily produced in lattice 
simulations. For this study two interpolating fields are 
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FIG. 3: Probability distributions for the amplitudes and 
masses in lattice units extracted from local correlators using 
Np = 2 Wilson fermions at pion mass 500 MeV on a lattice 
of spatial length 1.8 fm at /3 = 6.0 using Jjv. 



quality data. In addition, we observe that the state of 
lowest mass that is present in the mass spectrum of the 
correlator computed with Jn is absent when using J^. 
Instead the correlator with has a lowest state that 
does not show up when using Jj^ . The conjecture is that 
this state is the Roper. 



III. FORM FACTORS 

To extract information on hadron structure one needs 
to calculate coupling constants, such as the nucleon axial 
charge qa, the ttN and ttTVA coupling constants, form 
factors, moments of parton distribution functions and 
generalized form functions. In order to compute these 
quantities we need to calculate the relevant three-point 
functions, which, in addition to the forward propaga- 
tor needed for the calculation of the masses, require the 
evaluation of the sequential propagator. The three-point 
function, related to the matrix element of the operator 
O between hadron states \h' > and \h >, is given by 

(G'''<^'*(t2,ii;p',p;r)) - ^ exp(-zp' •X2)exp(+zq-xi) 





0.001 0.002 0.003 0.004 0.6 0.8 1 1.2 1.4 
AO MO 




FIG. 4: As in Fig. [3]but using the interpolating field 



considered: 



Jn{x) 



^abc 



(2) 



As can be seen from the histograms shown in Figs. [3] 
and|31 one clearly identifies the first excited state in the 
positive parity channel of the nucleon using rather low 
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where for O we consider the electromagnetic and axial 
currents. We use sequential inversions through the sink, 
which allows us to obtain the three-point function for any 
momentum transfer q and operator insertion but fixes the 
quantum numbers of the initial and final baryons. 



A. Nucleon Electromagnetic form factors 

The elastic nucleon electromagnetic form factors are 
fundamental quantities characterizing important features 
of neutron and proton structure that include their size, 
charge distribution and magnetization. An accurate de- 
termination of these quantities in lattice QCD is timely 
and important because of a new generation of precise 
experiments. The matrix element of interest is 
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where p{s) and p'{s') denote initial and final momenta 
(spins) and ttiat is the nucleon mass, Fi{0) — 1 for the 
proton and F2(0) measures the anomalous magnetic mo- 
ment. These form factors are connected to the electric, 
G_e(9^), and magnetic, Gm(9^), Sachs form factors by 
the relations 
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(5) 



To extract the nucleon matrix element from lattice mea- 
surements, we calculate, besides the three point function 
G^-'''^(t2, ii; P P; r), the nucleon two-point function, 
G^^{t, p), and look for a plateau in the large Euclidean 
time behavior of the ratio 



R{t2,ti;p ',p ;r;/i) = 



(G^^"^(f2,ti;p',p;r)) 
(G'^^(t2,p';r4)> 



1/2 



(G^^(t2-ti,p;r4)) (G^^(ti,p';r4)) (G^^(t2,p';r4)) 
_ (G^^(t2-ti,p';r4)) (G^^(ti,p;r4)> (G^^(t2,p;r4)> 

'^"*^>4''*^»'n(p',p;r;^) . (6) 

where 

(G'^'^(f,p;r)) =^e-"P-" r*^" {n| T J"(x,t)J''(O,0) |J7 ). 

X 

(7) 

We use the lattice conserved electromagnetic current, 
j^{x), symmetrized on site x and projection matrices for 
the Dirac indices 
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Throughout this work we use kinematics where the fi- 
nal nucleon state is produced at rest and therefore q = 
p' — p = —p. For the polarized matrix element one can 
construct an optimal linear combination for the nucleon 
sink, which in Euclidean time is given by 



-S'm(q; i) = ^ n(-q ■,Tk;ii = i) = 



C 
2m N 



{P2 - P3)Sl 



+ {P3-P1)S2,^ + [pi -P2)<53,. )Gm{Q^) (9) 



with = —q^. This construction provides the maxi- 
mal set of lattice measurements from which Gm{Q^) can 
be extracted requiring one sequential inversion. No such 
improvement is necessary for the unpolarized matrix el- 
ements given by 



n(0,-q;r4 ;/i 
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and 



n(0,-q;r4;M = 4) = C^^^±^ GEiQ') 



(10) 



(11) 



which yield Ge{Q^) with an additional sequential inver- 
sion. G = rl?-T7:-^ is a factor due to the normal- 
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FIG. 5: Nucleon isovector electric form factor using Nf = 2 
TMF. 
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FIG. 6: Nucleon isovector magnetic form factor using Np — 2 
TMF. 



Besides using an optimal nucleon source, the other im- 
portant ingredient in the extraction of the form factors is 
to take into account simultaneously in our analysis all the 
lattice momentum vectors that contribute to a given . 
This is done by solving the overcomplete set of equations 



P(q;/i) = i?(q;/i)-F(Q2) 



(12) 



the different sink types, F 



and D is an M x 



where /'(q; /i) are the lattice measurements of the ratio 
given in Eq. ^ having statistical errors Wk and using 

Ge 
Gm 

2 matrix which depends on kinematical factors with M 
being the number of current directions and momentum 
vectors contributing to a given Q^. We extract the form 
factors by minimizing 
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using the singular value decomposition of D. The anal- 
ysis described in this Section to extract Ge{Q^) and 
Gm{Q^) is also applied to the analysis of all form factors 
presented in this work. 

The J N ^ N transition contains isoscalar photon 
contributions. This means that disconnected loop di- 
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FIG. 7: Nucleon magnetic moment using Nf = 2 TMF (filled 
circles), TVf = 2 + 1 DWF [31] (crosses) and Nf ^ 2 Wilson 
fermions [32] (open triangles) . The physical value is shown by 
the star. 



agrams also contribute. These are generally difficult 
to evaluate accurately since the all-to-all quark propa- 
gator is required. In order to avoid disconnected dia- 
grams, we calculate the isovector form factors. Assuming 
SU{2) isospin symmetry, it follows that ( p|(|?27^u — 
iJ7^d)|p) - ( n|(|w7'^u - iJ7^d)|7i) = ( p\{uj''u - 
d'y^d)\p) and therefore by calculating the proton three- 
point function related to the matrix element of the 
right hand side of the above relation we obtain the 
isovector nucleon form factors G^{q^) — G^{q^) and 

The results for the isovector electric and magnetic form 
factors using Np = 2 twisted mass fermions are shown 
in Figs. [5] and [6] [1^. The lattice results on the electric 
form factor fall off slower as compared to a parametriza- 
tion of the experimental data [26| shown by the solid 
line, whereas the magnetic form factor is closer to ex- 
periment. This is consistent with recent high accuracy 
results obtained by the LHP Collaboration ^] and the 
RBC-UKQCD Collaboration [Hj. Fitting the magnetic 
form factor to a dipole form we extract Gm(0), which 
determines the anomalous magnetic moment. We show 
its dependence on the pion mass in Fig. [T] Using chiral 
effective theory with explicit nucleon and A degrees of 
freedom to one-loop order the isovector anomalous mag- 
netic moment j29| . the Dirac and Pauli radii can be ex- 
trapolated to the physical point [H, [2^, [s^l • There are 
three fit parameters for the magnetic moment and the 
best fit to the twisted mass data with the associated er- 
ror band is shown in Fig. [71 Multiplying the Pauli radius 
squared with the magnetic moment yields an expression 
with only one-parameter like the Dirac radius that can 
shift the curves but does not affect their slopes. As can 
be seen the physical magnetic moment is within the error 
band whereas for the radii results closer to the physical 
point are needed to check the predicted slope. 
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FIG. 8: The Dirac radius squared (top) and Pauli radius 
squared multiplied by the magnetic moment (bottom). The 
notation is the same as that of Fig. [T] 
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Nf ^2 + 1 DWF. 



B. Nucleon axial form factors 

The matrix element of the weak axial vector current 
between nucleon states can be written as 



1/2 



u{p,s') 



GAiq'h,l, + I^G,{q'] 



-u{P,s) (14) 



where the axial isovector current — ij){x)^ ^j_'y^^tp{x) . 
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FIG. 10: Top: The nucleon axial form factors Ga{Q^) for 
TVf = 2 TMF and A^F = 2 + 1 DWF. The dashed line is a 
dipole fit to the lattice data whereas the solid line to experi- 
ment. Bottom: Gp(Q^) for Nf ^ 2 TMF. The dashed line is 
predicted from Ga(Q^) and Eq. (|15p . The dotted lines is the 
best fit with the associated error band. 



fermions at a comparable value of the pion mass. The 
Q^-dependence of Ga{Q^) can be well described by a 
dipole Ansatz go/(Q^/w^ + 1)^ as shown by the dashed 
line. This is what is usually used to describe experimen- 
tal data for Ga{Q^) where a value of itia ~ 1-1 GeV 
is extracted for the axial mass. However the axial mass 
TTiA extracted from the lattice data is larger resulting in 
a slower fall off as compared to experiment shown by the 
solid line. Assuming pion pole dominance Gp{Q^) can be 
obtained in terms of Ga{Q^) as 



Am%/ml 
1 + Q'^/ml 



Ga{Q^)- 



(15) 



In Fig. [To] we show with the dashed line what pion pole 
dominance predicts if we use the fit determined from 
GAiQ^)- The error band shows the best fit to Gp{Q'^) if 
we instead fit the strength and mass of the monopole in 
Eq.[T5l 



IV. N TO A TRANSITION FORM FACTORS 

The determination of the N to A electromagnetic and 
axial transition form factors requires the evaluation of the 
three-point function {G^^^ {t2,ti]p ' ,p;r)) with a new 
set of inversions, where for O we consider the electro- 
magnetic and axial currents. The 7*A^A matrix element 
is given by 



Having computed the nucleon electromagnetic form fac- 
tors we can obtain the axial ones with no additional 
inversions [25l [3^ . The advantage here is that only the 
connected diagram contributes. In addition at zero mo- 
mentum transfer we obtain the nucleon axial charge gA, 
a quantity that is very accurately measured experimen- 
tally. We show in Fig. [9] results obtained usingA^i? = 2 
twisted mass [l^ and domain wall fermions j3l|. The 
leading one-loop chiral perturbation theory result for gA 
in the small scale expansion [34j can be used to extrapo- 
late lattice results to the physical point. Making a three- 
parameter fit to the twisted mass results we obtain the 
solid curve shown in Fig. [Altogether with the error band 
determined by allowing the fit parameters to vary within 
a increase by one unit from the minimum. Note that 
this error band does not include uncertainties in the fixed 
parameters. We obtain, at the physical point, a value 
with a large error that is lower than the experimental 
value by an amount slightly larger than one standard 
deviation. Results closer to the physical pion mass are 
needed to reduce the error due to the chiral extrapola- 
tion. 

The Q^-dependence of the nucleon axial form factors 
Ga{Q^) and Gp{Q'^) using Np = 2 twisted mass fermions 
is shown in Fig.lTOj Our results for Ga{Q^) are in agree- 
ment with those obtained using Np = 2-1-1 domain wall 
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The evaluation of the two subdominant electromagnetic 
form factors Ge2[Q^) and Gc2{Q^), which are of pri- 
mary interest as far as the question of deformation is 
concerned, require high accuracy. A lattice QCD calcu- 
lation accurate enough to exclude a zero value to one 
standard deviation would point to deformation in the 
nucleon/A system. This is particularly relevant given 
the fact that extraction of these form factors from exper- 
iment involves modeling and therefore a non-zero value 
from a first principles calculation even to one standard 
deviation is an important result. Optimized sinks are 
constructed to isolate the subdominant form factors [sH] 
along the same lines as discussed for the polarized nu- 
cleon matrix element. In experimental searches for de- 
formation, it is customary to quote the ratios of the elec- 
tric and Coulomb quadrupole amplitudes to the magnetic 
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FIG. 11: The EMR calculated using quenched (stars) and 
dynamical (filled circles) Wilson fermions and in the hybrid 
approach (squares). 
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FIG. 13; Ga/i(Q^) in the hybrid approach and using Nf = 
2 + 1 domain wall fermions. Quenched Wilson results are also 
included for comparison. 
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FIG. 12: The CMR calculated using Wilson fermions and in 
the hybrid approach. 



Results on these ratios obtained using Wilson fermions 
and the hybrid action are shown in Figs. fTTl and [T2l Lat- 
tice results at low are non-zero. The lattice values 
of CMR at small are less negative than experiment. 
As the pion mas decreases lattice results tend to become 
more negative approaching experiment. Therefore one 
anticipates that for even smaller pion masses the dis- 
crepancy between lattice and experiment will be reduced 
since pion cloud effects are expected to make CMR more 
negative as we approach the physical regime |36| . 



In Fig. [13] we compare results for the dipole form factor 
Gmi{Q^) obtained within the hybrid approach and using 
dynamical domain wall fermions at about the same mass. 
As can be seen there is very good agreement showing that 
results obtained within the non-unitary hybrid action are 
reliable. 

The invariant proton to A+ weak matrix element is 
expressed in terms of four transition form factors as 
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FIG. 14: The ratio of to A axial transition form factors 
C6iQ'^)/CiiQ'^)- The dotted line shows the pion pole dom- 
inance prediction for the hybrid case. The dashed and solid 
lines are fits to a monopole form for the hybrid and quenched 
results respectively. 
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, riA, 2-. , Cjj {q ) 
+<^5 (Q )5'Am H 2 — l^lf' 



up{p,s) (16) 



where C^{Q'^) and Cq{Q^), the dominant form fac- 
tors, can be related assuming pion pole dominance like 
Ga{Q^) and Gp{Q^) are related in the nucleon case. 



In Fig. M we plot the ratio C^(Q2)/C^(q2)^ rpj^g 
doted line shows the prediction assuming pion pole dom- 
inance after a dipole fit to the hybrid results on C^{Q^) 
is performed. As in the nucleon case the ratio does not 
fall off as rapidly and a fit to a monopole form to the 
same hybrid results gives the dashed line. 
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O _2 




□ quenched Wilson, m =410MeV 
hybrid, = 353 IVIeV 
^ dynamical Wilson, m_ = 384 MeV 



0.5 



in GeV^ 



1.5 



FIG. 15: The Q^-dependence of Ge2{Q'^)- The green (red) 
hne and error band show a dipole fit to the mixed action 
(quenched ) results. The value of Ge2, in units of e/{2m\), 
at = are —0.810 ± 291 for the quenched calculation 
-0.87 ± 67 for Nf = 2 Wilson case and 
hybrid calculation. 



2.06^2 35 for the 



V. A ELECTROMAGNETIC FORM FACTORS 
AND DENSITY DISTRIBUTION 

Since the A(1232) decays strongly, experiments to 
measure its form factors are harder than for the N to 
A transition and yield less precise results. The A form 
factors can be computed using lattice QCD more accu- 
rately than can be currently obtained from experiment. 
The decomposition for the on shell 7* A A matrix element 
is given by 

{A{pf,Sf)\j'^^\A{p,,Si))^Au,{pf,ss)0''^^Ur{p,,Si) 



2mA 



(17) 



where ai{q'^), 02(9^), ci(q^), and 02(9^) are known 
linear combinations of the electric charge form factor 
Geo{<1^), the magnetic dipole form factor Guiio^), the 
electric quadrupole form factor GE2{<f'), and the mag- 
netic octupole form factor Gm3{<1^)- An optimized sink 
is constructed that isolated the subdominant electric 
quadrupole form factor [3, [s^l ■ The results are shown 
in Fig. HSl 




FIG. 16: Quark transverse charge density in a A"*" polarized 
along the s-axis, with s± = -1-3/2. The light (dark) regions 
correspond with the largest (smallest) values of the density. 



frame. Fig. 1161 shows the transverse density p^^^ for a 
A+ with transverse spin s± — +3/2 calculated from the 
fit to the quenched Wilson lattice results for the A form 
factors (which has the smallest statistical errors of the 
three calculations). It is seen that the A+ quark charge 
density is elongated along the axis of the spin (prolate). 



VI. CONCLUSIONS 

Lattice QCD simulations are now being carried out in 
the chiral regime by a number of collaborations. We have 
shown that there is agreement among recent lattice re- 
sults using different fermion discretization schemes on the 
low lying baryon spectrum and the nucleon form factors. 
Recent results on the low lying hadron spectrum where 
lattice artifacts have been carefully examined are in per- 
fect agreement with experiment providing validation of 
the lattice approach and QCD itself [9]. Furthermore we 
have shown that lattice QCD provides a framework for 
the computation of quantities that can not be accurately 
measured in experiment such as the A form factors pro- 
viding valuable insight into the structure of such hadrons. 
We anticipate that other key hadronic quantities will be 
computed to sufficient accuracy and with lattice artifacts 
taken into account thereby providing direct comparison 
to experiment. 



The electric quadrupole form factor is particularly inter- 
esting because it can be related to the shape of a hadron. 
Just as the electric form factor for a spin 1/2 nucleon can 
be expressed precisely as the transverse Fourier trans- 
form of the transverse quark charge density in the infinite 
momentum frame [ssj , a proper field-theoretic interpre- 
tation of the shape of the A(1232) can be obtained by 
considering the quark transverse charge densities in this 



Acknowledgments 

I would like to thank my close collaborators T. Ko- 
rzec, G. Koutsou, J. W. Negele, C. N. Papanicolas, E. 
Stiliaris, A. Tsapalis and M. Vanderhaeghen for their 
valuable contributions. I would also like to thank my 



9 



collaborators in the ETMC for their important contri- 
butions in various aspects of the work presented here 
and for a most enjoyable and fruitful collaboration. This 
work is partly supported by the Cyprus Research Promo- 



tion Foimdation under contracts nENEK/ENIEX/0505- 
39, EPYAN/0506/08 and KY-rA/0907/11. 



[1] R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz [Alpha 

collaboration], JHEP 0108, 058 (2001). 
[2] R. Frezzotti and G. C. Rossi, JHEP 0408, 007 (2004). 
[3] C. Alexandrou et al. [ETM Collaboration], Phys. Rev. D 

78, 014509 (2008). 

[4] C. Alexandrou et al. [ETM Collaboration], PoS 
LAT2007, 087 (2007) arXiv:0710.1173 [hep-lat]]. 

[5] R. Frezzotti and G. Rossi, PoS LAT2007, 277 (2007). 

[6] A. Shindler, Phys. Rept. 461, 37 (2008). 

[7] M. Gockeler et al. [QCDSF Collaboration], PoS 
LAT2007, 041 (2007). 

[8] S. Aoki et al. [PACS-CS Collaboration], Phys. Rev. D 

79, 034503 (2009). 

[9] S. Durr et al, Science 322, 1224 (2008). 

[10] L. Del Debbio, L. Giusti, M. Luscher, R. Petronzio and 
N. Tantalo, J HEP 702, 056 (2007). 

[11] Y. Aoki et ffl/. larXiv :0903.4155 [hep-lat]. ^ 

[12] C. Aubin et al, Phys. Rev. D 70, 094505 (2004). 

[13] Ph. Hagler et al [LHPC Collaborations], Phys. Rev. D 
77, 094502 (2008). arXiv:0705.4295 [hep-lat]]. 

[14] C. Alexandrou et al, Nucl. Phys. A 825, 115 (2009). 

[15] P. Boyle [RBC Collaboration and UKQCD Collabora- 
tion], PoS LAT2007, 005 (2007); C. AUton et al [RBC- 
UKQCD Collaboration], Phys. Rev. D 78, 114509 (2008). 

[16] S. Aoki et al [JLQCD Collaboration], Phys. Rev. D 78, 
014508 (2008)^ 

[17] K. Jansen, la,rXiv:0810.5634l (hep-lat] . 

[18] C. Alexandrou, G. Koutsou, J. W. Negele and 
A. Tsapalis, Phys. Rev. D 74, 034508 (2006). 

[19] Drach et al, PoS LAT2008, 123 (2008) 

[Exrv:0905.2894' [hep-lat]]; C. Alexandrou et al, in 
preparation. 

[20] A. Walker-Loud, ef al, Phys Rev. D 79, 054502 (2009). 

[21] C. Michael, Nucl. Phys. B259, 58 (1985); M. Liischer 
and U. WolflF, Nucl. Phys. B339, 222 (1990). 

[22] S. Basak et al, Phys. Rev. D 76, 074504 (2007); B. G. 
Lasscock et al, Phys. Rev. D 76, 054510 (2007); C. Gat- 
tringer et al, Phys. Rev. D 78, 034501 (2008); C. Morn- 
ingstar, arXiv:0810.4448 [hep-lat]. 

[23] P. Lepage et al, Nucl. Phys. Proc. Suppl. 106 (2002) 12; 
C. Morningstar, Nucl. Phys. Proc. Suppl. 109A (2002) 
185; T. Draperet al, arXiv:0810.5512 [hep-lat]. 

[24] C. Alexandrou, E. Stiharis, C.N. Papanicolas, PoS 
LAT2008, arXiv:0810.3982 

[25] C. Alexandrou et al [ETMC], 'arXiv:0811.0724' [hep-lat]. 

[26] J. J. Kelly Phys. Rev. C 70, 068202 (2004). 

[27] J. D. Bratt et al [LHP Collaboration], PoS LAT2008, 
141 (2008) arXiv:0810.1933 [hep-lat]]. 

[28] S. Ohta and T. Yamazaki [for RBC and UKQCD Col- 
laborations and for RBC and UKQCD Collaborations], 
[axXrv :0810.0045 [hep-lat]. 

[29] T. R. Hemmert and W. Weise, Eur. Phys. J. A 15,487 
(2002). 

[30] M. Gockeler et. al, Phys. Rev. D 71, 034508 (2005). 
[31] T. Yamazaki et oL, arXiv:0904.2039 [hep-lat]. 



[32] C. Alexandrou, G. Koutsou, T. Leontiou, J. W. Negele 

and A. Tsapalis, Phys. Rev. D 76, 094511 (2007). 
[33] T. Yamazaki et al [RBC-I-UKQCD Collaboration], Phys. 

Rev. Lett. 100, 171602 (2008). 
[34] T. R. Hemmert, M. Procura and W. Weise, Phys. Rev. 

D 68, 075009 (2003). 
[35] C. Alexandrou et al, Phys. Rev. D 77, 085012 (2008). 
[36] V. Pascalutsa and M . Vanderhaeghen, Phys. Rev. Lett. 

95, 232001 (2005). 
[37] C. Alexandrou et al, Phys. Rev. D 79, 014507 (2009). 
[38] M. Burkardt, Phys. Rev. D 62, 071503 (2000) [Erratum- 

ibid. D 66, 119903 (2002)]. 



